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ON THE DESIGN OF AIRFOIIS IN WHICH THE TRANSITION 
OF THE BOUNDARY LAYER IS DEIAYED" 


By Itiro Tani 
INTRODUCTION - LAMINAR-FLOW АТКЕОТІ5 


1. In high speed flight conditions, the drag of an airfoil is 
almost exclusively due to skin friction. Therefore, if further reduc- 
tion in drag is desired, it is necessary to delay as much as possible 
the transition from laminar to turbulent flow in the boundary layer along 
the surface, thus decreasing the extent of the turbulent boundary layer 
which gives considerable skin friction. As the factors that may affect 
the transition, we will consider the stream turbulence, the surface 
roughness, the surface pressure distribution, and so on. In actual 
flight condltions, however, the effect of turbulence seems to be unex- 
pectedly small, so that, so far as smooth surfaces are concerned, there 
remains only the shape of the airfoil section in relation to pressure 
distribution as the most important factor affecting transition. Ме call 
a laminar-flow airfoil that airfoil in which the shape of the section 
is suitably designed so as to delay the transition of ihe boundary 
layer flow, : 


2. It is evident that the laminar separation of the boundary layer 
may cause the transition, as will be mentioned in the appendant part of 
the paper, paragraphs 35-40. We cannot expect, therefore, to maintain 
laminar flow beyond the separation point. Summarizing the results of 
flight experiments on airfoils hitherto made (refs. 8 to 12), we have 
the conclusion that the observed transition coincides approximately with 
the calculated laminar separation point at small Reynolds numbers, while 
it moves upstream toward the minimum pressure point as the Reynolds 
number increases. However, no example has ever yet been observed in 
which the transition moves ahead of the minimum pressure point. We 
therefore arrive at the supposition that the laminar-flow airfoil may be 
most simply realized by designing the airfoil in which the minimum pres- 
pure occurs well downstream. 





*"КубквА во no Sen'i o okuraseru Yokugata ni tuite.’ Report of the 
Aeronautical Research Institute, Tokyo Imperial University, No. 250 
(vol. 19, no. 1), Jan. 1943. 
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DESIGN ОЕ SYMMETRICAL ATRFOIIS IN WHICH THE MINIMUM 


PRESSURE OCCURS DOWNSTREAM 


3. Following Professor Moriya (ref. 13), we write the coordinate 
along the chord in the form x = = (1 + cog $), and assign x = O, 
Ë = n to the leading edge, апа x = 1, 6 = О to the trailing edge. 


Expressing the ordinate of the mean camber line by M = ) a, сов nË, 
О 


and the half-thickness measured normal to the chord by T = 3 ba sin ni E 
| the pressure distribution around the airfoil in the aedi bend 
potential flow is given by 


2 


cos a peint +b nan (1 - cos nÊ) sin nË 


Заг" 
1 
+ sin a - = (1 ~ cos D i) nan sin ab + ) my cos nt 


2 


г sin? 6 + е a пау, sin 08 + W nb, сов n$ 
1 1 


. + eB.) 








Тһе ordinates of the upper and lower surfaces are given by M + Т 
and M - Т, respectively. 
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where о is the angle of attack, p is the pressure acting on the air- 
foil surface, measured from the static pressure of the undisturbed stream, 
and а is the dynamic pressure of the undisturbed stream. We assume that 


the trailing edge is sharp, so that у nb, = 0. We limit the range of 
1 : 


the variables EË between O and л, and assign the upper and lower 
parts of the double sign for the upper and lower surfaces, respectively, 
Writing Ғ for the maximum value of M (the maximum camber) and e for 
the maximum value of 2T (the maximum thickness), we put 





= 1 - cos nŠ sin nÉ aM 
РА = -2 n РА = «Ὁ = = mm 
c NE sin Ё š: Es "sin b üx 
со sin nË Y cos nË ат 
d 1 віп Ё 6 E! n sin t | dx 


Тре lift coefficient is given by 


ο... пан COS а 


We consider first only the thickness of the airfoil (the camber of 
the center line will be considered in the next section, paragraphs 8 
to 11). Namely, we consider the. symmetrical airfoil section set at zero 
angle of attack, with a view to obtaining the minimum pressure well 
downstream, | 


l. We adopt as the typical example of the commonly used symmetrical 
airfoils the NACA symmetrical airfoil (ref. 13) 


T=e Ша - 0.6300х - 1.7580x7 + 1.42153 - 0.5075 } 


The maximum thickness is located at x = 0.3, the leading-edge radius is 
1.1e* Í| and the trailing-edge slope -(ат/ах) хє] is 1.1(е. Тһе pressure 
distribution for the саве е = 0,1 ів shown in figure 1. Тһе minimum 


pressure is located at х = 0.1, and the laminar separation point, 
determined by the approximate method due to the author (refs. 15 and 16), 
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at x = 0.61. If the transition point of the boundary layer would not 
move upstream beyond the laminar separation point, we might expect to 
maintain a laminar boundary layer for more than half the surface of the 
airfoil. Тһе flight experiments hitherto made, however, appear to give 
negative evidence for such a conjecture, 


5. Now, in order to shift the minimum pressure backward, it is 
required to shift the position of maximum thickness (x = m) backward. 


For designing such airfoils, we represent the shapes of parts before 
and after the maximum thickness by two algebraic expressions. For the 


forward half (0 S x S m) 


ZI D 


while for the rear half (m Ś x Š 1) 


T = 322 + dy(1 - х) + dp(1 - х)? + dz(1 - өз) 











where 
2 - зна {5πα - 1 
найган зас 2 
2m 
1.47 - 281(1 - m) d,(1 - m) - 0.98 
do = d4 = 3 
(1 - m)? (1 - m) 


and we assign arbitrary values for three parameters, m, В (= leading- 
edge radius 4. е2), and di (= trailing-edge slope +e). Although the 
method has the drawback that the two expressions give different values 


of а2т /ах2 at x = м, where dT/dx becomes zero, we nevertheless 
adopt it because we are in a position to vary the forward and rearward 
parts most simply and independently. 


6. First, we fix the forward half with m = 0.5 and h = 0.5, and 
vary the rear half by giving dy the values 1.7, 2.0, 2,5, and 3.0, 


respectively. Тһе shape of the section and the pressure distribution 
for е - 0.1 are shown in figure 2. We find from thie result that, ав 
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d4 increases, the minimum pressure point moves backward and the gradient 


of pressure rise following the minimum pressure steepens. We also find 
that the pressure distribution іп the neighborhood of the minimum pres- 
sure exhibits a wavy indentation when the value of dđ is too small or 


too large, and that there exists a certain value of dy for which the 


pressure distribution is flat and smooth. Such a value of dj is about 
2.5 in this case. We therefore fix the rear half with dj = 2.5, and 


vary the forward half by giving h the values 0.35, 0.50, 0.70, and 
1.05, respectively. The shape of the section and the pressure distri- . 
bution for е = 0.1 are shown in figure 3. From this comparison, we 
Pind that the negative pressure bump immediately behind the leading 
edge decreases as h decreases, and that the maximum permissible value 
of h is about Ο.Τ. 


The effect of thickness is shown in figure h, in which curves of 
pressure distribution are given for different values of e, 0.06, 0.10, 
and 0.14, but for a fixed set of parameters, m = 0.5, h = 0.5, and 
di = 2.5. It is seen that the characteristics of the pressure distri- 
bution do not materially change with thickness. There is, however, a 
slight change in the pressure distribution, the maximum permissible 
value for h slightly increasing as the thickness increases. 


To see the effect of the position of maximum thickness, we give m 
values ranging from 0.35 to 0.60, varying at the same time values of dy 
and h so that the pressure distribution becomes flat and smooth. The 
result of calculation is given in figure 5, which shows a considerable 
Change in the position of minimum pressure. The change is not purely 
due to the effect of m, but it is at any rate to be noticed that the 
value of m less than 0.4 is not sufficient for shifting backward the 
minimum pressure, while increasing the value of m beyond 0.5 is of no 
advantage, since the backward shift is then almost saturated, only the 
adverse pressure gradient being increased. 


7. From the results of calculation, we thus arrive at the conclusion 
that m must be between 0.4 and 0.5 and h must be less than 0.7 in order 
that the minimum pressure occurs well downstream. Smaller values of h 
are desirable, but, on the other hand, we should like to make Һ ав 
large as possible, because a large value of h will be advantageous in 
increasing the maximum lift coefficient and in preventing the inception 
of adverse pressure gradient when the angle of attack is slightly changed. 
Even if we give h the maximum permissible value 0.7, the leading-edge 
radius amounts to only 60 percent of that for the conventional NACA sym- 
metrical airfoil of the same thickness. Іп order to increase the leading- 
edge radius, it is required to increase the thickness, which in turn is 
accompanied by an increase in adverse pressure gradient following the 
minimum pressure. The adverse pressure gradient should be kept within a 
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certain limit, so that it becomes necessary to make a compromise between 
conflicting requirements. Thus, we are nọ longer in a position to 
require the farthest possible rearward location of the minimum pres- 
sure. We should also use a value of di which is somewhat smaller 


than that mentioned previously. 
Taking these requirements into account, we finally arrive at the 


design of a series of symmetrical airfoil sections, the parameters of 
which are given in the following table: 


Position of 
лы БЛ БЕГИ ESL 






Although section I is the most ideal for delaying the transition, in 
practice, its extraordinarily sharp nose and blunt tail are drawbacks. 

On the other hand, section N is too much compromised. Sections K or L 
seem to be suitable as laminar-flow airfoils for practical use. The 
ordinates of these six sections are given in table 1, while the auxiliary 
functions В, and Вс associated with the pressure distribution (see 


paragraph 3) are given in tables 2 and 3, respectively. The shapes of 
. the airfoil sections and the pressure distribution for е = 0.1 are 
shown in figure 6, 


DESIGN OF MEAN CAMBER LINE SUITABLE POR LAMINAR-FLOW AIRFOIIS 


8. A symmetrical airfoil set at zero angle of attack has no lift. 
In order to obtain lift, the center line of the symmetrical airfoil must 
be curved with a suitable camber. Since the effects of thickness and 
camber are nearly additive with regard to the pressure distribution, the 
mean camber line which maintains the nature of the pressure distribution 
of the symmetrical airfoil will be such that it shall give a uniform dis- 
tribution of pressure difference when the thickness is removed. Evidently, 
the center of pressure is then located at x = 0.5, so that such a camber 
line has the drawback that the travel of center of pressure is consider 
able. To reduce the travel of center of pressure, the uniformity of pres- 
sure difference should be satisfied only in the forward part of the chord. 
From the standpoint of designing the laminar-flow airfoil, however, it is 
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only required that the distribution of pressure difference is uniform 
from the leading edge to that point corresponding to the minimum pres- 
sure of the symmetrical airfoil. 


9. When the angle of attack a is small, the expression for pres- 
sure distribution given in paragraph 3 may be put into the form 


2 





Ë + eB) 2 2 "A = x a) - (ға, = eB. Ja 


2 
1 + (fA, £ eBo) 





ια [rd 
|| 
~ 
1 


Since the effect of the term (ға, + eBe} is very small, the quantity 


ІС = ТАС + 





is required to be constant in order that the camber line shall not change 
the nature of the pressure distribution of the symmetrical section. The 
range of constancy is at least up to the position of minimum pressure of 
the symmetrical section. Putting сов б = u = 2x - 1, and considering 
for simplicity the case when the minimum pressure is located at u =O, 
we prescribe that 


G = constant = Go for 0<х©0.5, -LŚUŚO 


G2 
| 


m 
Goll - ue), m20 for 0.5 &х 51, 02012201 


See figure 7. Moreover, since 


55 Fes 125 
Цэ” сов nb , cos 5 
l вір E gin 6 


a cannot be arbitrary, but must be so chosen that the right side of the 
equation does not become infinite at the leading edge, Š = sx. It is 
given by 4: 
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where (') denotes that only odd integers should be taken Рог п. This 
is the go-called ideal angle of attack due to Theodorsen (ref. 17). 
Using the assumption of the thin wing theory, we neglect the тегин eB, 


and (fA, + eB ë). We then have 


х/2 
Ст, = 2(1 + σι) 160 σι = 1 sin™tie at 


1 + (1 + ш) о, 


C2 
В 

ll 

' 


+(1 + m)(1l + og) 


со С 


E “L 
шы 25 nan + Zw ex 


where Cpo 18 the moment coefficient about the leading edge (positive 
when nose up) at Сү = O. Although m = О corresponds to the case of 
making G uniform up to the trailing edge, it seems to be impossible 


to realize a finite pressure difference at the trailing edge. Moreover, 
the quantity -Cyo/Cy, (which represents the degree of center of ргев- 


Sure travel) is as large as 0.25 in this case. If ш»0, G vanishes 
at the trailing edge, and -Сио/Ст, decreases ав m increases, tending 


to О as m approaches о, ш = о corresponds to the case when G = 0 
in the rear half of the chord. Increasing the value of m, however, 
steepens the pressure gradient, so the value of m from 3 to 5 seems to 
be adequate, 
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Now, since 


4 CI, Сү, 
fG sin Ë =2 Y na, cos nË - [2 na, + — | сов Ё + — 
1 1 2х 2x 


the slope of the camber line having the prescribed distribution of 6 
is given by 


dM $^ sin nb) 1: ΠΤ у ағ 
ах 1. sin 6 n JO 1L : z сов § - сов ΗΝ 








ыз ie - log Sie ey” δε δε 
ох 1 + On 


penc а У ш 


u- vV 1 


* 


The ordinate of the camber line is obtained by the integration 


co 
Y na, may be determined by the condition that М = 0 at x= 1. We 
1 


call Ру the camber line thus determined. The equations for camber 
lines for m= 0, 1, 3, 5, and о, namely Do, Di; 03, 55, and D,, 
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are given below, their important characteristics being summarized in the 
table. It is to be noted that f is the maximum value of М, є is а 
the absolute value of the zero-lift angle, and є .and а are measured. 









in radians. 2 
НЕСЕТ 
010 | 
. «0380 
.0609 
Dg: M vj ocz: 2 (b ua) PEPE: 10 Q3) 
f 2 log £ 5 | 
| n 
0: gp M = (5 + u) lg 2 - — (1 + u) log (1 + u) + u? log jul - 


= (1 - u) (2 + u) log (1 - u) «Ξ (1 - uê) 


Lex 1 35 | 
3* GF 51 (51 + 19u) log 2 51 (1 + u) log (l + u) + 


1. 1, 
Я u3(35 - 2145 + 5и ) log u | - 
= ας u)*(16 + 29u + 2002 + 543) log (1 - u) + 


— (1 - u£) (176 - Blu - 172:2 + 30u3 + бом) 
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D5: 25 M = εἷς (949 + Ата) iog 2 - ŚŚ (1 + u log (1 + u) + | 
sip «(1155 - 1386:2 + 99048 - 38506 + 6348) log ja] - 
Xs 0 = u) (256 + Виза + 121842 + 93803 + 378ut + 


6302) log (1 - — (i-a - 28535u - 66088u2 
зы?) log (1 - u) в" 42) (35072 кете ын + 


—— 


3168003 + 687924" - 1743092 - 3612006 + 3780uf + 15608) 





D: M {а + u) log 2 - (1 + u) log (1 + u) + и log lal | 
f log 3 - | 


Shapes of these camber lines are shown in figure 7, the ordinates 
of them are given in table №, and the auxiliary functions Ас and Ag 
(see paragraph 3) and the pressure difference distribution G are given 
in tables 5, 6, and 7, respectively. | 


10. The calculation made previously is only approximate, neglecting 
the thickness. І% is therefore desirable to check the result by actually 
calculating the pressure distribution for the specified angle of attack 
taking both camber and thickness into accomt, Ав an example, we con- 
struct an airfoil by applying the thickness form К with е = 0,15 τοτ- 
mal to the chord around the camber line Ds with f = 0.02 ( the 
resulting airfoil is designated as D=K - 2015). We calculate the pres- 


sure distribution by the formula of paragraph 3 for the optimum design 
condition а = 0.99? and Ст = 0.245. The result is shown in figure 8. 


The nature of the pressure distribution remains similar to that of the 
symmetrical airfoil, so we may consider that the approximate determina- 
tion neglecting thickness gives results sufficiently accurate for practi- 
cal purposes. | | Ён 


ll. In designing the camber line Од, we have assumed for simplicity 


that the pressure difference G is constant for u < 0. This corresponds 
to the case when the symmetrical airfoil has its minimum pressure in the _ 
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neighborhood of u = О. Therefore, the camber line Dy is adequate to +. - 
be combined with the symmetrical section 2 or К. 17, however, the 
symmetrical section is adopted in which the minimum pressure is located 
further upstream, it is not only not necessary to maintain G constant 
up to u= О, but also of disadvantage because it makes it difficult to 
reduce the value of  -Cmo f. 

To reduce the range over which С should be maintained constant, 
we may proceed in the following way. Assuming for instance that G 


should be constant from u= -l to u=- =, and using a new variable 


uj = Е (1 + 3u), ме prescribe that 


G = Go for -1 © us - $ ~ = 5 450 
a Gą(1 2)3 1-2 < < < 
= Goll - ul Рог 5 weny 0 š u š 1 


The calculation may be performed similarly to the case of Dy. The . ' 
resulting camber line is designated ав F3. The camber line lying in 
the middle between F3 and Ds is also designed, and designated ав Ej. 


Their Important characteristics are given in the following table together 
with those of Dz. Other numerical data for these camber lines are given 


іп tables # to 7. К 





Since the camber lines ZE) and Fa enable us to maintain G con- Dr 
stant up to the point x - 0,52 and x - 0.33, respectively, they are 
adequate to be combined with the symmetrical sections Г and M, respec- 
tively. The pressure distribution is shown in figure 10 for the airfoil 
obtained by applying the thickness form M with е = 0,15 around the 
camber line Ёз with f = 0.02. The optimum design condition corresponds 


to а = 1.100 and Ст = 0.252. 
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EXPERIMENTS ON LAMINAR-FLOW AIRFOILS 


12. In order to ascertain whether it is possible to prevent the 
forward movement of the boundary layer transition by shifting the mini- 
mum pressure on the airfoil surface, we have to perform experiments in 
a low turbulence wind tunnel or on the actual airplane in flight. When 
the Reynolds number is not too large, however, we can still use a con- 
ventional wind tunnel in which the stream turbulence is relatively small. 
So we made at first comparative measurements on two symmetrical airfoils, 
NACA 0010 and L.B. 24 in the 1.5 m wind tunnel of the Aeronautical 
Research Institute. І.В. 24 is a laminar-flow airfoil of 10 percent 
thickness, already shown in figure 3. The theoretical pressure distri- 
bution is also given in figure 11. The minimum pressure is located at 
х = 0.64, and the laminar separation at х = 0.77. The wind tunnel was 
of the lowest turbulence level available for the author, the critical ` 


Reynolds number of & sphere being 3.66 x 102 and the transition Reynolds 


number of a flat plate 1.05 X 103 (see paragraph 28). In order to raise 
the Reynolds number as high as possible, unusually large models were 
used. They were made of laminated mahogany, of highly polished surface, 
of 0.8 m span, of 1.2 m chord, and fitted with end plates 1.3 m X 0.6 m. 
Since the model was large compared with the size of the tunnel and the 
end plates were not sufficiently large, the results for a given airfoil 
may not correspond even approximately with those for the same airfoil in 
an undisturbed two-dimensional flow. Our object, however, was merely 

to ascertain the relation between pressure distribution and transition, 
and it seemed reasonable to expect that the relation will not be seriously 
affected by limitations in the conditions of the experiments. As a matter 
of fact, marked difference was found in the calculated and measured dis- 
tributions of pressure, the latter of which was measured along the median 
section of the model with a static tube of 1 mm diameter (fig. 12).2 

This discrepancy, however, is immaterial, since our object was merely to 
compare the two airfoils, both of which are affected quite gimilarly by 
experimental limitations. 


13. The angle of attack of the model was zero , and the wind speed 
was varied from 6 to hO m/s. Тһе local drag of the median section was 
determined from wake measurements, that were made in the section 11 cm 
behind the trailing edge. Measurements of static and total pressures in 
the wake were made, respectively, with a static tube of 2.5 mm external 
diameter and a pitot tube with a Plattened mouth of 0.65 mm external 
depth and 2.6 mm width. The profile drag coefficient Cpy was obtained 





“The measured values are those for a Reynolds number of about 2 X 106, 
The distribution of pressure changes but little with the Reynolds number, 


Ф 
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from the measured pressures by Jones‘ formula (ref. 18). Figure 13 3 
presents Cp, plotted against Reynolds number R referred to chord 


length. For a lower range of R, the drag of L.B. oh is higher than 

that of NACA 0010, the reason probably being that a turbulent boundary 
layer associated with the laminar separation 1в established at a higher 
Reynolds number for the former airfoil than for the latter. For a higher 
range of R, however, the condition is reversed, L.B. 24 giving a drag 


less than half that of NACA 0010 for R higher than 2 x 106. This is 
probably due to the fact that the transition may occur much later for 
Т.В. 24 than it does for МАСА 0010, as also observed from the compari- 


son of wake conditions for the two airfoils (fig. 14) 4 


14, In order to verify the aforementioned supposition, a pitot tube 
with a flattened mouth of external depth 0,9 mm and width 2, 1 mm мав 
placed in contact with the airfoil surface, and the wind speed, and con- 
sequently the Reynolds number R, were determined at which the indicated 
total pressure 0% divided by the dynamical pressure а of the undis- 
turbed stream begins to rise suddenly. The results are shown in fig- 
ure 15. From this figure, the dependence of the transition point on 
Reynolds number as shown in figure 16 is obtained. At the same Reynolds 
number, the transition occurs much farther from the leading edge for 
L.B. ой than for NACA 0010. Even at the ша Reynolds number reached, 
L.B. 2} has a transition ав far back as = 0.80. This is somewhat 
beyond the laminar separation point, x = ο 1, which 18 calculated from 
the theoretical pressure distribution. However, this is not contradictory, 
because the actual pressure distribution differs from the theoretical one 
in a manner to delay the transition (fig, 12). 


15. With further increase in Reynolds number, the transition шау 
move toward the leading edge, but it seems improbable that the transition 
moves forward beyond the minimum pressure. 1% is highly desirable to 
check this point also by wind tunnel experiments, but all the wind tunnels 
now available to the author are of no use for making measurements at 
sufficiently high Reynolds numbers, because the transition is prematurely 





JIn this figure, the curves L and Т represent the drag of a 
flat plate when the boundary layer is entirely laminar &nd entirely 
turbulent, respectively. Тһе curves NV and NF represent the drag 
Of airfoil NACA 0009 measured in the NACA Variable-Density Wind Tunnel 
and NACA Full-Scale Wind Tunnel, respectively, 


Ат this figure G and p are the total and ЖЕТЕ pressures іп 
the wake, respectively, and Gg ів the total pressure outside the wake, 


all being. measured from the static pressure of the undisturbed stream. 
у is the distance across the wake, and. % is the chord of the model. 
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induced by the turbulence of the stream (see paragraphs 27 to 29). It 
seems urgent to build a special wind tunnel of low turbulence level. 
For the present, however, it is simplest.to rely upon experiments in 
actual flight. Such a hope of the author was fortunately realized by 
the specially planned flight experiment, which was performed at the 
Navy Aeronautical Technical Arsenal (ref. 19). 


16. The airplane used for the experiment was a biplane; two portions 
of the lower wing, each of 1.1 m span, were covered with the airfoil to | 
be tested. The test portions were of chord 2.4 m, made of Japanese Hinoki, 
highly polished, and fitted with a partition fence of small height at 
both ends. Two test portions were placed symmetrically, pressure distri- 
bution and wake measurements being performed on the starboard portion, 
while the boundary layer was observed on the port portion. Тһе airfoil 
section was not one of the most sppropriate design now considered, because 
lt was required to put it on the original section of the sirplane, and, 
moreover, to determine the section before completion of the final design 
calculation. 1% has the following characteristics: 


Mean camber line: M = 0.0667х(1 - х) (1 - x + х2), f = 0.0125 


Thickness distribution: е = 0.12, m = 0.45, h = 0.56, dy = 1.60 


The camber line is similar to Dg of paragraph 9, but there exists a 


slight lack of uniformity of G in the neighborhood of leading and 
trailing edges. The thickness distribution is similar to L of para- 
graph T, but the trailing edge slope is somewhat larger than  L. 


lT. Results of flight experiments are summarized in figure 17 
and 18. In figure 17, the section lift coefficient Ст, obtained by 
integrating the pressure distribution curve, is shown by a broken line 
plotted against the Reynolds number R referred to the flight speed 
and chord length, and Ст is again shown by a solid line plotted against 


the profile drag coefficient Сто determined from the wake measurements, 


In figure 18, the measured pressure distribution is shown in comparison 
with the theoretical one (two-dimensional potential flow) having the 
same value of Cy. The transition points estimated from the change in 
boundary layer velocity profiles are also marked. Generally speaking, 
the measured pressure distribution agrees fairly well with the theo- 
retical one, although a slight difference appears when Ст, becomes 
large. An adverse pressure gradient is found on the lower surface when 
Ст, is small, thus resulting in the transition point being observed 
unexpectedly far forward. Such a discrepancy in pressure distribution ` 
as compared with the theoretical one seems to be probably due to the 
fact that the span of the test portion was not sufficiently large. As 
a result, the profile drag coefficient [Do has the minimum value 0.0042 
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at about. Ст, = 0.26, which is larger than the value Ст = 0.18 theo- 


retically estimated on the assumption that the transition occurs far 
back on both upper and lower surfaces, Therefore, the observed value 
of Cpg, although much smaller than that of the conventional airfoils, 


seems to be still somewhat large when compared with the optimum case. 

At any rate, however, no transition was found to occur upstream of the 
minimum pressure. It is important to note that such experimental evi- 
dence was obtained on an airfoil section in which the minimum pressure 

is located further downstream than on the conventional one. This finding 
will give valuable data to establish a basis for design of the laminar- 
flow airfoils. 


ESTIMATION OF THE DRAG OF IAMINAR-FIOW AIRFOIIS 


18. As mentioned previously, the results of flight experiments 
seem to support the basis for the design of laminar-flow airfoils, 
namely, the possibility of maintaining the boundary layer laminar at 
least up to the minimum pressure point. It is interésting, therefore, 
to estimate the drag of laminar-flow airfoils by assuming a laminar 
boundary layer from the leading edge to the minimum pressure point and 
a turbulent boundary layer downstream to the trailing edge. 


For the laminar boundary layer, the momentum thickness is given by 
o O.hhy ΓΡ h 
Ө = ua. ds 
0 


uy? 





with a sufficient approximation (ref. 16), where w ів the velocity 


outside the boundary layer and в ів the distance measured along the 
airfoil surface from the forward stagnation point. Writing t for the 
chord length and V for the velocity of the undisturbed flow (velocity 
of flight), and putting 


we have the nondimensional expression > 


~ 
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where the subscript a refers to the point of minimum pressure. 
Applying then the solution due to Buri (ref. 20) for a turbulent 
boundary layer assumed to extend from the minimum pressure point to 
the trailing edge, we have the result | 


PH iT κα - 0.016eR71/ ы |. = ut ac 


ба 


where the subscript b refers to the trailing edge. The numerical 
values originally given by Buri are slightly modified во as to agree 
with measurements when applied to the flat plate. 


According to Squire and Young (ref. 21), the profile drag coeffi- 
cient is given by 


CDQ = = (9 + өз) 03" 


Where the subscripts u and 1 refer to the upper and lower surfaces, 
respectively. Тһе exponent 3.2 of ІЛ, has been obtained by assuming 
the ratio of displacement and momentum thicknesses equal to 1.4. But the 
ratio seems to exceed 1.4 near the trailing edge, so we replace 3.2 by 


3.4 with a view to improving the accuracy and at the same time to 
simplifying the algebra. Since eet ee с 


Ө Ө 17/4 Е c 4/5 
Du a (Ge) T^ + 0.0162R "I Pag ο) 
| а 
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we have 


Cp, = 0.074872 -— А a р (1,2-3/8 " ар 


where 


5/8 | 

σ С: 

ΝΟ t= [ас 
σ 


Β 


ТР the velocity distribution uj = UV ів calculated by assuming the 


potential flow of an ideal fluid, it is desirable to modify the distri- 
bution to take account of the effect of separation near the trailing 
edge. We tentatively modified the distribution of U such that the 


Buri parameter T = (e fm) (виз уаз) (316 ~) + at the trailing edge for 


the case when the boundary layer is assumed turbulent from the leading 
edge, namely  - 


2 
| r _ 0.0081 (=) f^ piss 
6 do / Jo 


shall not become smaller than -0.06. In almost all the cases, values 
of Up thus modified are found in the range between 0.95 and 1.00. 


l9. Applying this method of calculation, the profile drag coeffi- 


cient Сто is estimated first for a series of -symmetrical airfoils set 


at zero angle of attack. The series consists of the six symmetrical 
airfoils, I, J, K, L, M, N, as given in paragraph 7 and the NACA con- 


ventional airfoil. Values of στο at ЕВ = 2 х 10Í for three different 


thicknesses (maximum thickness in terms of chord e =.0.10, 0.15, 0.20) 

are shown in figure 19 plotted against the position of minimum pressure, 

Сто seems to decrease almost linearly as -the minimum pressure is shifted . 
backward, the most ideal airfoil I giving a value about half of that of 

the NACA conventional airfoil. If it is desired to realize a profile 

drag of two-thirds of the conventional airfoil, it will be required to 

use the symmetrical airfoil L with the maximum thickness located at 

45 percent chord from the leading edge. | 
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Effect of camber is relatively small. ІТ, for instance, the 
center line of the symmetrical airfoil К with е = 0.15 is curved 
into the camber line Ds with f = 0.02 (see paragraph 9), the esti- 


mated increase in profile drag at the optimum angle of attack is only 
0.0001. For f= 0.04, it is 0.0003. 


Finally, we compare the laminar-flow airfoil with the most exten- 
sively used airfoil, NACA 23012, for which the leading-edge radius ів 
0.0158, and the optimum lift coefficient corresponding to the minimum 
profile drag coefficient is about 0.15. If we consider the symmetrical 
airfoil section К (h = 0.56) combined with the camber line Ds, it is 


necessary to use the thickness е - 0.15 in order to obtain the same 
magnitude of leading-edge radius, and the camber f = 0.012 in order to 
realize the optimum Lift coefficient 0.15.2 Therefore we construct an 
airfoil by applying the symmetrical form К with е = 0.15 normal to 
the chord around the camber line D5 with f=0. 012. We call it 


Dek - 1215. The angle of attack corresponding to Ст = 0.15 ig 1.569 


for МАСА 23012 and 0.60" for DsK - 1215. The pressure distribution for 
that condition is shown in figure 20. | SE 


—— — 


We then estimate CDO for the two airfoils by the method explained 


previously. The results are shown by broken lines in figure 21. In 
order to check the results, measured values taken from various sources 
for the two airfoils and similar airfoils are also plotted in the ваше 
Pigure by different marks. The mark o refers to the value obtained 

by Flight experiments on a smooth surface, and e refers to that 
obtained by wind tunnel experiments where the stream turbulence has no 
effect on transition. The mark + refers to the flight experiment on 

a rough surface, while X refers to the wind tunnel experiment where 
the stream turbulence causes the transition to occur prematurely. There- 
fore, only o and e are adequate for our present purpose. Drawing 
curves through these points and extrapolating to higher Reynolds numbers, 
we find that the result agrees fairly well with the estimated values. 
Therefore, we may consider that the method of estimating [Do is suffi- 


ciently accurate at the Reynolds numbers corresponding to actual flight 





ZThe calculation developed in paragraphs 8 to 11 refers to the poten- 
tial flow of an ideal fluid, so that it gives the slope of lift curve Ἢ 
С | 
— = Ол. In real fluids, however, the slope of lift curve amounts to 
a 
only 80 to 90 percent of the theoretical value. If we take this effect 
into account, we have to increase the necessary amount of Ї by 10 to 
20 percent in order to realize the given lift coefficient. However, 
such a slight change in the value of Ғ will scarcely affect the ев%1- 
mation of στο: | 
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conditions. Comparison of two airfoils, laminar-flow and conventional, 
also suggests the possibility of 40 percent reduction in profile SIUE 
by using а fairly practical laminar-flow airfoil. 


CONSIDERATION OF THE AIRFOIL WITH UNIFORM DISTRIBUTION OF PRESSURE 


20. The fact that, so fsr as flight experiments with smooth wings 
are concerned, the boundary layer transition occurs only in the region 
of rising pressure, not only warrants the principle of designing the 
laminar-flow airfoil by shifting the minimum pressure backward, but 
also suggests the possibility of delaying the transition by using an 
airfoil with uniform distribution of pressure. Therefore, in para- 
‘graphs 21 to 23, the shape of such а symmetrical airfoil is determined 
by a method similar to that used for designing the camber line of laminar- 
flow airfoils, and the airfoil was examined by wind tunnel experimenta. 
In paragraphs 24 to 26, a calculation is made to inquire about the method 
of sucking away the boundary layer over the region of rising pressure in 
Buch a way that the boundary layer velocity profile shall remain the same 
as that for the point of minimum pressure. 


21. Consider the symmetrical airfoil set at zero angle of attack. 


According to the formula of paragraph 3, the pressure distribution is 
given by 


| 2 
со віп пЁ 
1 + 2 у nb, 
[ gin 6 


с 
1+ |2 nb, ER 
i 


E-1- 
езі 








sin & 


where p is the pressure acting on the airfoil surface, measured from 
the static pressure of the undisturbed stream, а is the dynamic pres- 


sure of the undisturbed stream, x= 5 (1 + cos &) is the coordinate 


along the chord, and the half-thickness of the airfoil is expressed in 
the form 


co 


Рад. b, sin në 
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If the thickness is sufficiently small, the square of 


= cos në 
аг == «Ὁ > nbp ee сыган наз 
dx 1 sin $ 


may be neglected, so that the condition of uniform distribution of pres- 
sure is satisfied by putting all the coefficients b,, other than bj, 


equal to zero, namely, by an elliptic section. Іп order to take the 
thickness into account approximately, we substitute the value of dT/àx 
for the elliptic section into the denominator of the expression for p/q. 
Then, writing e for the maximum thickness in terms of the chord, we 
get 


2 


co 


sin Ë + 20 nb, sin nt 
I 





віп2ё + efcos £ 
hence 


2 > nb, sin nË = - sin Š + BJsin2E + e*cos2$ 
L 


1/2. 
where В is the constant value of (1 - p/q) / . We have therefore 
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т ©з ОР: 4 2245: шшш m 1 
3 in Ë E 
sin 6 = = - z / 2 > nb, sin nŠ ` — a = 
О 1 сов Ё ~ сов 5! 
x t t 
= JJ sin 8 - ТЕЛ + е©сов®®' b _ min $' 8 
0 cos Ё ~ сов Εξ) 
= ú - = 2ku'cog le + = 
55 (1 + u) (1 - ku + е\1 - x&v?) 
1 - kuć log > 
(1 - u) (1 + k2u + ер 5 keue) 
where 


u = cos Ё = Px - 1 k = γι - eż— 


Upon integrating we get | | 2 | 23 


Т = f x cosmo ДРЕА. | Е u 


> EE ы (1 u u) (1 - Bu + el - 22) du 
ex Jat V 1 - u2 (1 - а)(1 + kêu + еуі - xeu?) 


The integral is evaluated by a numerical method, and the value of the 
constant В determined from the condition that T = że when и = 0. 
The numerical results for three values of е are given in table 8, Т 
and aT/dx being expressed in terms of those for the elliptic section. 
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The shapes of the airfoils are shown in figure 22. The shape resembles 
an elliptic section, although it is somewhat fuller at the ends. It will 
also be seen that the constant B, as shown in the following table, is 


slightly smaller than 1 + e, the maximum value of (1 - р/а)? for the 


elliptic section. Since the values of (dT/dx) + (-eu | l - 42) are not far 
different from 1, it seems to be sufficiently accurate to substitute 

the value of 4Т/4х for the elliptic section into the denominator 

of 1 - p/q. 





22. The uniform distribution of pressure requires, however, an 
infinite pressure gradient at both leading and trailing edges. Іп order 
to see to what degree such a sharp pressure gradient may be realized in 
actual fluids, measurements were made on & model of the airfoil section 
with uniform distribution of pressure with е - 0,1 (we call it U.P. 0010) 
in the 1.5 m wind tunnel of the Aeronautical Research Institute, The 
model was made of laminated mahogany, of 0.8 m span, of 0.8 m chord, and 
fitted with end plates 1.3 m X 0.6 m. Measurements of pressure distribu- 
tion, wake traverse and boundary layer transition were similar to those 
already mentioned in paragraphs 12 to 1}. 


The pressure distribution along the chord is shown in figure 23 
for three values of R, the Reynolds number referred to chord length. 
Ihe observed value is somewhat high compared to the theoretical 
value = = -0.203, the discrepancy probably being due to the excessive 
size of the model in proportion to that of the wind tunnel, At any rate, 
however, the pressure distribution is nearly uniform. Тһе lack of uni- 
formity exists at both edges due to the impossibility of realizing the 
infinite pressure gradient. The boundary layer separates near the Ἢ 
trailing edge, but the effect of separation becomes small as the Reynolds 
number increases. This scale effect seems to be of the same nature as 
that responsible for the sudden drop in sphere drag; the boundary layer 
separates in a laminar state when the Reynolds number is low, while it 
becomes turbulent before separation when the Reynolds number is high, 
thus being able to proceed against a larger pressure gradient. This is 
also seen from the measurements in the wake, where the indentation of the 
curve of total pressure distribution is shallow and wide for low Reynolds 
numbers, while it becomes deep and narrow as the Reynolds number increases. 
As a result, the profile drag coefficient στο decreases considerably 
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as the Reynolds number increases, as shown in figure 24. Transition to 
turbulence was found very near to the trailing edge, occurring downstream 
of х= 0,9 in the range of measurements. Measured values of the drag 
of the model when a piano wire of 0.5 mm diameter was placed at x = 0.8 
and x = 0.9, respectively, are also plotted in the same figure. The. 
drop in drag occurs at a lower value of the Reynolds number when the 
surface is roughened by the wire. 


23. The profile drag coefficient of the airfoil U.P. 0010 is shown 
below in comparison with other symmetrical airfoils at R = 2.2 X 106; 


МАСА 0010 Сро = 0.0064 (fig. 13) 
Г.В. 2} | 0.0032 (fig. 13) 
U.P. 0010 0.0059 
О.Р. 0010, a wire at х = 0.9 0.001 


All the airfoils are of 10 percent thickness. NACA 0010 is a conventional 
airfoil, and L.B. 24 is a laminar-flow airfoil with far back minimum 
pressure. The drag of U.P. 0010 is between that of these two airfoils, 
the drag when a wire is placed being nearly the mean of the two. This 
result seems to be interesting in that the drag of an airfoil with a 
blunt tail is smaller than commonly considered. 


The airfoil with uniform distribution of pressure will also probably 
be favorable when used at high subsonic speeds. Even if the shock wave 
occurs at high subsonic speeds, the increase in drag will remain small 
when the boundary layer does not separate. This expectation was really 
verified by the experiment due to Kawada and Kawamura (ref. 22), the 
drag of the airfoil U.P. 0010 being smaller at high Mach numbers as 
compared with other airfoils. 


21, From the fact that the boundary layer transition occurs only in 
the region of rising pressure, we may also expect suction of the boundary 
layer to delay transition. For example, if the boundary layer is sucked 
into a slot, there is a well-known sink effect (ref. 23) which relieves 
the adverse pressure gradient somewhat upstream of the slot. We may con- 
sider an alternative possibility. That is, we assume that the boundary 
layer is laminar in the region of falling pressure, and that it remains 
laminar also in the region of rising pressure provided that the boundary 
layer profile is the same as that at the minimum pressure point. We 
then ask what suction arrangement must be applied in order to realize 
such a condition. 


25. We denote by s the coordinate measured along the surface, 
y perpendicular to the surface, 6 the boundary layer thickness, u 
the velocity in the boundary layer, uj the velocity outside the boundary 


4R 
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d: 
layer, == = -р (du; /as) the pressure gradient, and то = -H (du/dy) y=0 


the skin friction at the surface, Assuming the surface (y = 0) is made 
porous, through which the fluid is sucked with the velocity c, we have 
the equation of continuity | 


“-°-/ чау (а) 


а 
UW - — u = To + Š — 
PUH р : ду О is - (b) 


where w ів the velocity of fluid entering the boundary layer through 
y = 5. The equation of motion reduces to 


а 2 
-C >) = чу — + ν =) (c) 
7 у= Š у? y=0 


for y = 0. 


Now, the velocity profile in the boundary layer may be approximated 
by the Pohlhausen polynomial 


u = u 28-203, у (а) 


when neither pressure gradient пог suction exists. For this profile we 
have | 


21 
ΤΟ = Сы 75 | | | (е) 
Wh 
pu; = —— To (f) 
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To simplify the calculation, we assume that the expressions (d), (e), 
and (f) still hold when both pressure gradient and suction exist. Then 
we have from (b 


as _ 630 v _ 10h 5 21 (а) 


m 52 | т 2) 1260 p? аа 
u — = [11 + — u ds (h) 
1 № 1 8-80 37 t 


208 | 
where m = 367 snd во ів the initial position of suction, which is 


the minimum pressure point in the present case. We have also from (a) 


„ 1841 ὅπ κ (4) 
3670 ds 


If we substitute (d) and (1) into both sides of the equation (c), however, 
the left and right sides become 1.003; (du; /ds ) and u; (duy /ds), геврес- 


tively. This contradiction is evidently due to the crude assumption of 
using (f) in spite of the presence of pressure gradient and suction, but 
we may overlook the error because it is small. 


26. We apply the calculation to the symmetrical laminar-flow air- 
foil of 10 percent maximum thickness, І.В, 24, set at zero angle of 
attack. The velocity distribution 11/7 calculated from the potential 


flow of ideal fluids is used, 6 the maximum velocity (minimum pressure) 
being located at 64 percent ОР the chord from the leading edge (во - 0, 65ο). 


Applying a distributed suction downstream of the minimum pressure point 
во as to maintain the velocity profile in the boundary layer the same as 
at this point, we have the boundary layer thickness 9 and the required 





ότι is assumed that the velocity distribution is not affected by 
the suction. Theoretically u] should be О at the trailing edge, but 
the distribution was somewhat modified so as to give uq = 0.85V there. 
The effects of these assumptions appear to be too small Бо affect the 
result materially. 
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suction velocity c as shown in figure 25. It is to be noted that s 
is measured along the surface from the leading edge, and R is the 

Reynolds number based on the chord length t and the velocity of the 
undisturbed stream V. Integrating the area under the curve of c, we 


have the total amount of suction 1.3, VR per unit span of the two sur- 
faces, If we assume span = 35 m, t = 5m, V = 200 m/s, V = 0.15 cm^/s, 
the total amount of suction amounts to 5.6 m3/s, which will require an 
exit area of only 0.028 m^ when discharged with the velocity equal to V. 
integration of то gives the drag coefficient στο = 0.0005. This value 
may be compared with Cp, = 0.0003 for the flat plate with laminar 
boundary layer, στο = 0.0044 for the flat plate with turbulent boundary 
layer, and Сро = 0.0025 for І.В. 2% without suction. If the thickness 


of the airfoil is doubled (20 percent chord), then the amount of suction 
Will be nearly doubled; the drag is however almost unchanged. 


It should be noticed again that the calculation is based on the 
assumption that no transition occurs if the velocity profile in the 
boundary layer maintains the form at the minimum pressure. 1% is the 
purpose of the calculation to show that extraordinarily low profile drag 
may be expected with a relatively small amount of suction under such a 
condition. 


PREMATURE TRANSITION OF BOUNDARY LAYER ~ EFFECT OF STREAM TURBULENCE 


27. Although the transition of the boundary layer occurs only down- 
stream of the minimum pressure, so far as flight experiments on smooth 
airfoil surfaces are concerned, there are many examples of wind tunnel 
experiments in which the transition moves upstream of the minimum pres- 
sure. This seems to be due to the premature transition caused by the 
turbulence in the wind tunnel stream, For example, the transition on 
the airfoil L.B. 24 was found only downstream of the minimum pressure 
in the range of Reynolds numbers covered by the author's wind tunnel 


experiments (the Reynolds number based on chord length up to 3 X 106; 

see paragraphs 12 to 14); as a result very low values of the profile drag 
coefficient Сро were observed, The same airfoil, however, when tested 
with a larger model of 2 m chord in the 2.5 m wind tunnel of Kawasaki 
Aircraft Company, Gihu, gave the result as shown in figure 26, in which 


στο increases considerably when the Reynolds number exceeds 3 X 106 


(refs. 24 апа 25). There is reason to believe that the increase in drag 
is due to the effect of stream turbulence. The boundary layer observa- 7 
tion at the Kawasaki wind tumnel shows that the transition is found at 


50 percent chord (x = 0.5) for the Reynolds number 6 X 106 and moves 
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further forward as the Reynolds number increases. The boundary layer 


x< 
velocity profile observed at transition has a form factor => = 2.6 


(8% апа Ө are the displacement and momentum thicknesses of the 
boundary layer), which is very near to the value for the case of zero 
pressure gradient, This result seems to suggest that the transition is 
at least not correlated with the laminar separation (ref. 1). 


28, In order to verify this conjecture, it is desirable to show 
that the transition in the boundary layer along a flat plate occurs 
under the same condition, because that transition may be considered to 
be independent of the laminar separation. Unfortunately, however, no 
flat plate was measured in the Kawasaki wind tunnel. Therefore we pro- 
ceed in a somewhat indirect way. We assume that the degree of stream 
turbulence is represented by the conventional critical Reynolds number 
of the sphere, Rc, and the condition of transition due to turbulence 


u8 
represented by the local Reynolds number, Rg = =, at transition on a 


flat plate, where uj is the velocity outside the boundary layer and 

6 is the momentum thickness of the boundary layer. It is generally 
accepted that the turbulence in the wind tunnel stream will give a 
fluctuation of pressure gradient, as a result of which an instantaneous 
and intermittent separation will occur. Such an instantaneous and inter- 
mittent separation, however, does not necessarily lead to the transition 
into turbulence; for the transition really to occur, it seems probably 
necessary that the Reynolds number Rg which represents the ratio of 


inertia pu, to viscous stress рушт /θ exceed & certain critical value. 


lt 1s also expected that the critical value depends on the degree of tur- 
bulence; it must increase ав Rọ increases. This is really shown by 

the experimental data hitherto published, which are given in the following 
table and also by white circles (o) in figure 27. The available data are 
scanty, especially because the experiment on a flat plate is very diffi- 
cult. It was necessary for the author to perform a new experiment (ref. 30) 
with a view to adding one point in the range of high Ёс. 


26 


0.21 x 103 | National Bureau of Standards 
“το х 103 | National Bureau of Standards 


42 X 103 | N.P.L. Compressed Air Tunnel 
3,66 x 102 |1,05 x 103 | Aero. Res. Inst. 1.5 m Tunnel 
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29. Now, we calculate the value of Кө at transition of the а1г- 
foil 1.28. 21 tested in the Kawasaki wind tunnel, and correlate it with 
the critical Reynolds | number Rc of that tunnel, We analyze similarly 
the other available data, and summarize the result in the following 
table. The values of the form factor 87/0, not shown, were all found 
in the range from 2.1 to 2.7. We then plot the data also in figure 27 
by black circles (6). The black circles are seen to define a single 
curve together with the white circles already mentioned. This result 
seems to support the supposition that the transition under consideration 
is mainly caused by the stream turbulence, but not correlated with the 
laminar separation. 


1.85 х 10210.11 x 103|Symmetrical &irfoillN.P.L. 7 ft 
2.10 x 107| .56 x 103|Airship model MIT. 75 ft 


3.50 X 102| .95 x 103|Airfoil N-22 NACA Full-Scale 
3.65 x 102|1.08 x 103|Airfoil Т.В. 24 Kawasaki 2.5 m 






















In reference 12 (the third line in the preceding table), the same 
airfoil was examined both by the full-scale wind tunnel and by the flight 
tests. We calculate the form factor 8%/Ө from these tests and plot the 
values against s/t in figure 28, where t is the chord length and s 
is the length measured along the surface from the leading edge. The 
value of 5*/0 at transition is 2.6 in the wind tunnel, while it reaches 
as high as 3.1 and drops sharply in the flight test. The minimum ргев- 


sure is located at = = 0.18, and the laminar separation calculated from 


the measured dlstribution of pressure at = = 0.36. This example 18 very 


interesting because the cause of transition is quite different in the two 
cases (namely, it is due to the stream turbulence in the wind tunnel, 
while it is related to laminar separation in the flight test), although 
the positions of transition are almost the same, 


PREMATURE TRANSITION OF BOUNDARY LAYER - EFFECT OF SURFACE ROUGHNESS 


30. Up to this point, we have only considered the case when the sur- 
face of the airfoil is smooth. ТФ the surface is rough, however, there 
is a possibility that the transition may also be caused prematurely by 
surface roughness. So, it is important in practice to estimate the 
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approximate order of magnitude of the permissible roughness in the lami- 
nar boundary layer. Nothing has been known concerning this problem, 
except a шеге conjecture or fragmentary data. Schiller (ref, 36) sug- 
gested that a local separation occurs and hence leads to transition 

when the Reynolds number kuy/V exceeds a certain critical value Rrit’ 


where К is the height of projection and uy 18 the velocity at the 


top of projection. The exact value of Roerig 18 not known, but it is 
not likely to differ much from the critical value of the Reynolds number, 
above which vortices are shed from the obstacle of the same shape as that 
of the projection placed in a uniform stream, The experimental result 
due to Wieselsberger (ref. 37) shows that such a critical Reynolds number 
is roughly 50 for a circular cylinder. Assuming that the height of pro- 
jection k is small, and that the presence of the projection in no way 
alters the character of the flow, we have the shearing stress at the 


surface то = pV(ux/k). Using the so-called friction velocity Ух = (το/ο. 


instead of we have then kv,/V = |ku,/v., The permissible ro ess 
Ч x/ Ὡς Nm 


is therefore given by kvy/V = ҮН сүү, or, with Буу = 50, kvy/V= T. 
On the other hand, according to Nikuradse's experiments on roughened 


pipes (ref. 38), the critical Reynolds number is τ e 4, above which 


the roughness projections disturb the laminar sublayer of the turbulent 
boundary layer and hence increase the pressure drop. It appears there- 
fore that the permissible roughness is smaller in the turbulent boundary 
layer as compared with the laminar boundary layer. This is confined by 
a British flight experiment (ref. 39) on the airfoil section of 10 feet 


chord (Reynolds number 1.8 x 107), because the effect of camouflage paint 
of 0,001 inch thickness increased the drag by about 6 percent without 
moving the transition forward. At any rate, however, such a estimate is 
nothing but mere conjecture. With a view to making the estimate more 
definite, we performed wind tunnel experiments, although of small scale, 
of quantitative character (paragraphs 31 to 32}. 


31. A polished aluminum plate, 80 cm long, 60 сш wide, and 3 ши thick, 
was held horizontally in the 1.5 m wind tunnel of the Aeronautical Research 
Institute. So that the flow at entry would not be disturbed, the leading 
edge of the plate was rounded, and the plate slightly tilted so that the 
forward stagnation point was on the same surface as that where the observa- 
tion was made. The tilting, however, was so slight that the static pres- 
pure was observed to be practically uniform along the plate. The plate 
was roughened by a wire, which was stretched across the flow, in contact 
with the plate. The diameters k of the wire were 0.25, 0.4, and 0.7 m, 
respectively, and the distances x of the wire from the leading edge were 
15, 30, 45, and 60 cm, respectively. When the wind speed V was low, the 
boundary layer was laminar all along the plate, but from a certain speed 
upward, the transition to turbulent flow was observed at that point where 
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the wire was placed. Transition was detected by а sudden change in the 
value of the total pressure G* in terms of the dynamic pressure q of 
the undisturbed stream, G* being indicated by a small pitot tube with 
a flattened mouth of 1 mm external width and 0.3 mm width, which was 
placed in contact with the plate at a point 70 cm behind the leading 
edge. A sample record of measurements for k = 0.4 mm is shown in 


Pigure 29, 


When a flat plate is placed along a uniform stream of velocity V, 
the Blasius solution (ref, 40) of the laminar boundary layer equation 
gives 


= 0. στο E) A 


for a point of distance x from the leading edge. Writing K. for the 
critical value of kvx/v, the permissible height of projection К 1в 


given by 


κ _ z(VxN3/* 
0.576 = = (FP) 


т 


We determine У from the wind speed corresponding to the kink of the 
curve as shown in figure 29, plot 0.576(k/x) in a logarithmic scale 
against Үх/У, and draw a straight line of the slope -3/№ through the 
points (fig. 30). We thus obtain К = 13, which is far greater than 
the value K = 7 estimated previously. 


32, Similar measurements were also performed on an airfoil section 
Т.В. 21. Тһе model was of 0.8 m span, of 1.2 m chord, fitted with end 
plates 1.3 m X 0.6 m, and set at zero angle of attack in the same wind 
tunnel. Wires of various diameters (k = 0.25, 0.4, 0.7 mm) were attached 
parallel to the span, in contact with the surface, "at lO percent of the 
chord from the leading edge (x - 0.1). Transition was detected by the 
sudden change in total pressure as indicated by a pitot tube with a 
flattened mouth of 2.7 mm external width and 0.9 mm depth, which was 
placed in contact with the surface at СО percent of the chord from the 
a yt edge. Results of measurements are shown in figure. 31, where 


= 13, and t is the chord length. 


The friction velocity may be generally expressed in the form 
ув-1/% 
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where A is a function of s/t (в is the length measured along the 
surface from the forward stagnation point). We can calculate A by 
applying either the Pohlhausen approximate solution (ref. 41) or the 
simplified method due to the author (refs. 15 and 16). Тһе permissible 
height of projection is then determined by 


k -3/h 
A Š = КЕ 


Applying the Pohlhausen method to the theoretically calculated 
distribution of pressure, we get the values of function A, as shown 
in figure 32. Since A is 1.23 at the position of wire (x = 0.1), 
the value 1.23(k/t) is plotted in a logarithmic scale against R in 
figure 33, R being the Reynolds number corresponding to the kink of 
the curve as given in figure 31. It will be seen that although the 
measured points ere оп a straight line of the slope -3/h, they give 
K = 15, which is somewhat higher than the value found for the flat plate. 


33. Now we apply the preceding result to the fragmentary data 
hitherto known in order to check the adequacy of the estimate. First, 
we examine the results of wind tunnel experiment on a symmetrical 
laminar-flow airfoil L.B, ати on which various projections are attached 
at 3 percent of the chord from the leading edge (x = 0.03). The model 
was of 0.8 m span, of 1.2 m chord, and set at zero angle of attack. 

The profile drag was measured by the method similar to that for L.B. 24 
(see paragraphs 12 to 13). The results are shown in figure 34, from 
which we find that the rubber tape of 0.07 mm thickness gives no effect 
over the range of Reynolds numbers R covered by the experiment, while 
the piano wire of 0.5 mm diameter gives a completely turbulent friction. 
The effect of the wire of 0.25 mm diameter begins to appear at 


R= 1.3 X 106. Inserting the values’ К = 0.25 m, t=1.2m and 
А = 1.95 in the formula А £ = KR"3/*, we get К = 15.5, which is in 


good agreement with the result in paragraph 32. The value А = 1.95 
was read from figure 32, since the leading edge portion of L.B. 27 
almost coincides with that of І.В. 2k. Profile drag coefficients of 
both airfoils are also the same over the range of Reynolds numbers 
examined. 


34, In connection with the determination of boundary layer transi- 
tion on airfoils in the NACA full-scale wind tunnel (ref. 42), an aux- 
iliary measurement has been reported, in which the effect was examined' 





11.8. 27:has а maximum thickness of 10 percent of the chord at 
60 percent of the chord from the leading edge. See figure 5 of paragraph 6. 


NACA ТМ 1351 33 


of rubber tapes attached at 5 percent of the chord (x = 0.05) from the 
leading edge of the model. The airfoil section was МАСА 0012, the chord | 


was 72 inches, and the Reynolds number R мав 1,18 х 106, No effect 
was found when the tape was 0.003 inch thick, some effect began to 
appear when it was 0.006 inch thick, and the transition moved right to 
that position where the tape of 0.009 inch thickness was attached. 
Assuming A = 1.6 and K = 15, we estimate from the preceding formula 
the value 0.007 inch for the permissible thickness, which seems to agree 
well with the observation. 


If we further assume that the value А = 1.6 is also applicable 
to the case of the British flight experiment mentioned in paragraph 30, 
we find 0.004 inch for the critical height for transition with 


$ = 10 feet, R = 1.8 X 107, and К = 15. On the other hand, we esti- 


ку 
mate the permissible limit in turbulent boundary layer by E = В, 
"which шау be written in the form | КҮЗІ 


by the relation v = = сеш“, where uj is the local wind speed and 


Ср is the coefficient of local skin friction. If ме assume uj = 1.2V ᾿ =: 
and ср = 0.003 (which is equal to the coefficient of mean skin friction 


` for a flat plate at R = 1.8 х 101), we obtain 0,0006 inch for the per- ` 
missible roughness thickness. Since the thickness of the camouflage 
paint is reported to be about 0.001 inch, it may be concluded that the 
paint increases only the friction in turbulent boundary layer, without 
affecting, however, the transition to turbulent flow. This is in good 
agreement with the experimental results. 


34. | | NACA TM 1351 
APPENDIX 
TRANSITION CAUSED BY LAMINAR SEPARATION 


35. As is well-known, the phenomenon of sudden decrease in drag of 


a sphere at a certain value of the Reynolds number В = < (V is the 


speed of undisturbed stream and d is the diameter of sphere) is 
explained by supposing that the boundary layer separates while it is 
laminar when R is low, but it separates after transition to turbulence 
when R is high, thus resulting in diminishing the so-called dead water 
region. Probably the transformation from laminar separation to turbulent 
separation may proceed as follows: 


When the laminar boundary layer separates from the surface, the 
detached layer remains also laminar at first, but it is so unstable that 
it becomes turbulent at a short distance. This transition from laminar 
to turbulent flow is considered to occur when the local Reynolds number 
based on the width of the detached layer and the velocity outside the 
layer exceeds a certain value, so that the transition moves upstream 
toward the separation point as R increases. When the transition 
approaches sufficiently near the separation point, it becomes possible 
for the detached layer to come back again.to the downstream surface, 
because the turbulence produced will drive the flow forward. The layer 
reattaches to the surface as a turbulent layer, and accordingly the drag 
coefficient begins to decrease. The distance between the separation and 
the first turbulent boundary layer decreases as R increases, and finally 
the fully developed turbulent boundary layer commences just downstream 
of the separation point. The drag coefficient then ceases to decrease. 


36. Now, in order that the separated layer reattach to the surface, 
046 
it seems necessary for the local Reynolds number Rg = -4- at separation 


to exceed a certain critical value, where uq is the velocity outside 


the boundary layer, and 0 is the momentum thickness of the boundary 
layer. This may be explained as follows: According to the laminar 
boundary layer theory, the separation occurs when the quantity 


o du . wu 
paz R ŚP ϱ > py — 
Ү ав ds 8 


exceeds a certain value, suggesting that the pressure rise (dp/ds)9 
becomes too large in proportion to the shearing stress at the sur- 
face pv(uj/0). Assuming analogically that the separated layer leaves 
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the surface when the pressure rise becomes too large in proportion to 
the momentum puj”, we then find that Rg at the separation point must 


exceed a certain critical value in order that the кратер layer 
reattach to the surface. 


37. In order to determine the critical value of Вр, we consider 


in detail the condition where the coefficient of sphere drag begins to 
decrease. This condition corresponds to the point В of the curve оғ 
figure 35, which represents an idealized variation of the drag coeffi- 


cient p|" drag + (ov?/2) ο] or the pressure difference coeffi- 


cient Δδρ/α]- difference of pressures at the forward stagnation point 


and the point corresponding to the central angle 157.59 + (<v?/2) | with 


the Reynolds number R. Within the range AB, the pressure distribution 
around the sphere is approximately independant of R; the typical example 
may be found from the experiments due to Fage (ref, 13), Fortunately, 
the boumdary layer calculation has also been performed Рог that distri- 
bution of pressure by Tomotika and Imai (ref, 44), so that the local 
Reynolds number Rg is given by 


Bg = 0.40VR 


at the separation point. Although the calculation has originally been 


made for a particular Reynolds number, В = 1.57 X 102, the preceding 
relation may be applied for any value of R in the range of AB. Putting 
the value of В at В, and writing 


= š R Эн 1 R OZ . 
E о.о В 


ме һауе Ro orit as the critical value of Rg above which the separated 


laminar layer reattaches to the surface. Conventionally the Reynolds 
number Вс corresponding to Ср = 0.3 (or Ap/q = 1.22) has been used, 
instead of Rp, for representing the degree of stream turbulence, but it 
is not so difficult to estimate the value of Rp from the measured curve 
of Cp (or Др/а) against R. For example, we have from the experiments 
of towing spheres in the free atmosphere (ref. 35) 

= 240. 


Rp = 3.6 X 10? В = 3.85 x 105 


“бла 
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Since these values refer to the case of very low turbulence, it will be 
seen that HO ss pe 240 represents the highest possible value, It is 


R 
also to be noted that E 0.94 in this example and that almost the 


same value has been obtained by the author's experiments on spheres of 
various diameters (ref. 29). 


38. If the stream turbulence is not low so that Rọ is less than 


3.85 X 102, then Ra orit will be less than 240, Assuming the ratio Rg/Rc 
to be constant, we can estimate the corresponding value by 





On the other hand, we can also estimate the value of Rocr directly 


from the boundary layer measurements. The results of the analysis Рог а 
sphere as well as circular and elliptic cylinders are summarized in the 
following table, where Ro. „ję 18 the critical value Rg estimated 


from Вс by the nrecadine- formula, nO aap is the value of Rg observed 


at the separation point when the boundary layer really separates while 
it is laminar, and “буцаа: is the value of Rg observed at the calcu- 


lated laminar separation point when the boundary layer separates after 
transition. The fact that Ro... lies between 2277 and "O eau 


seems to suggest the adequacy of the preceding consideration. 


Гэ үх eo 


190 160 220 
150 140 225 С 
200 160 400 16, 47 










Sphere 
Circular cylinder 
Elliptical cylinder 2.7 x 102 















39. Ме пом proceed to apply our result to interpreting the effect ог 
Reynolds number on maximum lift of airfoils. For the angle of attack 
near the stall, the flow separates shortly downstream of the leading edge 
while the boundary layer is laminar. Jf the flow fails to reattach to 
the surface as a turbulent layer, the maximum lift coefficient Ошак of 
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the airfoil will be almost independent of the Reynolds number R, At 
the separation point, similar to the case of a sphere, the relation of 
the form - 


Во = KVR 


holds, where R = ΤΕ is the Reynolds number referred to chord length 5,8 


and К is a constant depending on the shape of airfoil and the value of 


μας" If В is low so that Rg is less than the critical value ἘΠ 


then Bez will be independent of R. Assuming the same value of 


PT? for the sphere as for the airfoil, we obtain 
Ke 


for the Reynolds number above which Ої ах begins to increase with the 


Reynolds number. Therefore, the ratio of the Reynolds number corresponding 
to a certain value о? πες of an airfoil and the critical Reynolds 
number of sphere in the same stream, Rp or Rc, becomes independent of 


the stream turbulence. Denoting the values for a reference tunnel with 
asterisk, we have 


R* = 
R - 


aja 


which in turn means that the ratio of Reynolds numbers corresponding to 
a certain value of πες is equal to the ratio of critical Reynolds 


numbers of a sphere. This is useful for comparing the values of C sss 


obtained in two different wind tumnels. Considering the reference condi- 
tion to be the free flight in the atmosphere, we find | 


ж 
R* = R x το 
Ro 





τε is to be noted that R is referred to t, while Rp and Rc 
are referred to d. 
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ав the free flight Reynolds number which will give the same value о? 
Ch ав that observed in a wind tunnel. This is just what is called 


the effective Reynolds number. Strictly speaking, such an argument as 
mentioned before should apply only to the Reynolds numbers near the 
critical value, but there are many experimental evidences showing the 
usefulness of the concept of effective Reynolds number for most practi- 
cal purposes, as far as the commonly used airfoils and range of Reynolds 
numbers of both wind tumnel and free flight are concerned. 


10. Finally, we consider a more quantitative example to show the 
adequacy of the preceding argument. In figure 36, “гах for various 


МАСА symmetrical airfoils are plotted against the effective Reynolds 
number R*, the experimental data being taken from the results of the 
NACA variable-dengity wind tunnel (ref; 48). As already mentioned, up ` 
to a certain value of R*, Claas is almost independent of R*, This 


corregponds to the condition in which the laminar separation just behind 
the leading edge fails to reattach to the surface, resulting in a con- 
siderable dead water region above the airfoil surface. The value of 

OT 88 is approximately 0.9, irrespective of the thickness; it is almost 


equal to the value for a flat plate of vanishing thickness. Theory of 
discontinous flow, when applied to the flat plate, seems to give a lift 
coefficient close to 0.9 (ref. 49). We idealize, therefore, the experi- 
mental curve ав shown by dotted lines in figure 36. Then, the point 
where the dotted line meets the line C * = 0.9 will be considered to 


correspond to Ro rit = 050. In order to determine this point, we calcu- 


late the value of Ro / at the laminar separation point for a lift 
coefficient Сү, = 0.9, and the value of R which gives Rg = 240. We 
first calculate the pressure distribution by the formula of paragraph 3 
for the two-dimensional potential flow around the airfoil section. 
Although the formula may be applied to any arbitrary airfoil section, 
we have determined the pressure distribution only for the airfoil 


= 0.28 Te lx (1 - x)(5 - hx) 


in order to simplify the calculation, because no great exactitude is 
required in the present problem. x is the coordinate along the chord, 
х= О and х = 1 corresponding to the leading and trailing edges, 
respectively, Т is the half-thickness, and e is the maximum thick- 
ness in terms of chord length. The airfoil represented by the preceding 
expression coincides with sufficient accuracy with the true NACA symmetri- 


cal airfoil except near the trailing edge. “The value of Ro / {В at the 


laminar separation point was then determined for the calculated pressure 
distribution by applying the approximate method due to the author (refs. 15 
and 16). 
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Тһе critical values R. ate thus calculated are shown by a solid 


line in figure 37, while the corresponding values taken from figure 36 
are shown by white circles. The agreement is fairly good, and especially 
satisfactory when the thickness of the airfoil is small. Іп general, the 
thickness of the boundary layer near the trailing edge increases as Cr, 


increases. If Ст is further increased, however, a laminar separation 


suddenly occurs near the leading edge when the thickness is small, while 
the trailing-edge turbulent separation moves a considerable extent for- 
ward before the leading-edge laminar separation occurs when the thickness 
is moderate. Therefore, the assumption of the analysis is more satis- 
factorily realized in the case of small thickness, thus bringing the 
calculated and observed values in close agreement. 


In conclusion, the author wishes to acknowledge his indebtedness 
for the assistance given by Messrs. С. Noda, Б. Mituisi, I. Shinra, 
5. Asaka, R. Hama, and К, Takeda. 


Translation by Itiro Tani 
University of Tokyo 
Tokyo, Japan 
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TABLE 1.- ORDINATES FOR SYMMETRICAL SECTION. VALUES OF T/e 
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TABLE 2.- AUXILIARY FUNCTION В, ASSOCIATED WITH THE 


PRESSURE DISTRIBUTION OF SYMMETRICAL SECTIONS 
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TABLE 3.- AUXILIARY FUNCTION В. ASSOCIATED WITH THE 


PRESSURE DISTRIBUTION OF SYMMETRICAL SECTIONS 
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TABLE 1.- ORDINATES FOR MEAN CAMBER LINES. 


VALUES OF M/f 
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TABLE 5.- AUXILIARY FUNCTION Ас ASSOCIATED WITH THE 


PRESSURE DISTRIBUTION OF MEAN CAMBER LINES 
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TABLE 6.- AUXILIARY FUNCTION Ag ASSOCIATED WITH THE 


PRESSURE DISTRIBUTION OF MEAN CAMBER LINES 
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TABLE 7,- PRESSURE DIFFERENCE DISTRIBUTION G FOR MEAN CAMBER LINES 
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TABLE 8.- AIRFOIL WITH UNIFORM DISTRIBUTION OF PRESSURE 
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